Coherence in complex networks of oscillators 
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We study fully synchronized (coherent) states in complex networks of chaotic oscillators, reviewing the an- 
alytical approach of determining the stability conditions for synchronizability and comparing them with nu- 
merical criteria. As an example, we present detailed results for networks of chaotic logistic maps having three 
different scale-free topologies: random scale-free topology, deterministic pseudo-fractal scale-free network and 
Apollonian network. For random scale-free topology we find that the lower boundary of the synchronizability 
region scales approximately as where k is the outgoing connectivity and /j, depends on the local nonlin- 
earity. For deterministic scale-free networks coherence is observed only when the coupling is heterogeneous, 
namely when it is proportional to some power of the neighbor connectivity. In all cases, stability conditions are 
determined from the eigenvalue spectrum of the Laplacian matrix and agree well with numerical results based 
on histograms of coherent states in parameter space. Additionally, we show that almost everywhere in the syn- 
chronizability region the basin of attraction of coherent states fills the entire phase space, and that the transition 
to coherence is of first-order. 



I. THE INTERPLAY BETWEEN DYNAMICS AND 
TOPOLOGY 



The structure and dynamics underlying complex networks 
have been widely investigated, providing insight for many sys- 
tems where they arise naturally GIBS]- Complex networks 
appear in a wide variety of fields ranging from lasers @|, gran- 
ular media quantum transport |7], colloidal suspensions 
0, electrical circuits 0], and time series analysis ficjl. to 
heart rhythms Hill, epidemics fl2l [l3 |j . protein folding I14tl . 
and locomotion Il5ll among others [JQyll- 

From the mathematical point of view, a network is a graph, 
composed by nodes or vertices and by their connections or 
edges |2|. Sometimes, each node is characterized by some dy- 
namical state (amplitude), which evolves according to some 
local contribution and to the interaction with the neighbor- 
hood. In other words, the complexity of the system underly- 
ing the network may be introduced either in the way nodes are 
interconnected (topology) or in the way nodes evolve in time 
(dynamics). 

When studying network dynamics one frequently assumes 
a regular topology (lattice) where each node evolves ac- 
cording to some more or less com plic ated dynamics, typi- 
cally fixed points fl6ll . limit cycles Il7ll or chaotic attractors 
One main goal of this approach is to study the 
so-called spatio-temporal chaos which appears in many dif- 
ferent spatially extended systems out from equilibrium, such 
as hydrodynamical flows, chemical reactions and biological 
systems l'20l l2lll . Two main topics in this context concern the 
study of mechanisms underlying pattern formation and pat- 
tern selection 1 2(1 l2lt l22l l23t uM and also the study of chaotic 
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synchronization behavior llri \2i 

Spatially extended systems are fundamentally modeled by 
(i) sets of coupled differential equations (20ll with nonlinear 
terms, where both time and amplitude are continuous, (ii) cel- 
lular automata ll22ll . where both time and amplitude are dis- 
crete or (iii) coupled map lattices l2lll . where time is discrete 
as in cellular automata, but the space of states is continuous. 
In all these models the underlying networks have connections 
whose range assumes all values from 1 (nearest neighbors) up 
to some maximum range, in particular the size of the system 
(global coupling regime). In other words, neglecting bound- 
ary conditions, these network systems assume translational 
symmetry and therefore the underlying network is called a 
regular network. 

To study more complicated network structures, one usually 
neglects node dynamics and all complexity is introduced by 
the network topology, i.e. by the way nodes are connected to 
each other. This may be done in three different ways Q|: by 
randomly connecting the nodes (random networks I26l 12711 ). 
by considering some random long-range connections in a reg- 
ular network with some small range of couplings (small-world 
networks I2H l29ll ). or by considering the introduction of 
new nodes which are connected to the previous ones follow- 
ing some rule of preferential attachment (scale-free networks 
IBOII '). For all these cases there is no translational symmetry 
and no typical range connection: connections do not have 
a 'spatial' interpretation. Therefore, one uses some general 
topological quantities to characterize each particular network, 
namely the average path length (£), i.e. the average minimum 
number of connections linking two nodes, the clustering coef- 
ficient C, defined as the average fraction of neighbors which 
are connected to each other, and the distribution of connec- 
tions P(k), representing the number of nodes having k con- 
nections. Table|l]shows the values of these three quantities for 
all three topologies. 

Random networks were introduced by Erdos and Renyi in 
the late fifties ll26ll to study organizing principles underlying 
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Random Small-world Scale-free 

(£) hi N/ In (pN) N for small p In N/ In In N 
In iV for large p 

C k/N Co(l-p) 3 ~N~ 3/4 

P(fc) e^P/kl eT l k k /k\ 2m 2 /k 3 

TABLE I: Characterizing complex topologies with the topological 
quantities: average path length {£), clustering coefficient C and dis- 
tributions P(k) of connections k. Here N is the total number of 
nodes, p is the probability for two nodes to be connected, k is the 
average number of connections per node, Co is the clustering coeffi- 
cient of the regular network from which the small-world network is 
constructed, and m is the number of initial connections of each new 
node in a scale-free network. 

some real networks izl . In random networks one defines 
some probability p(N), function of the total number N of 
nodes, which determines the probability for any two nodes 
to be connected in a total of N nodes. Consequently, the con- 
nections are typically long-range connections having a com- 
pletely irregular structure. One main goal in studying ran- 
dom networks is to determine the critical probability p c (N), 
beyond which some specific property starts to be very likely 
to be observed, e.g. the critical probability marking a tran- 
sition to percolation 13111 . One important feature of random 
networks, which also appears in real networks, is their small 
average path length (£), i.e. the average distance between any 
two nodes increases slowly with the system size. However, 
unlike random networks, real networks also have large cluster 
coefficients C. 

Small-world networks were introduced recently by Watts 
and Strogatz 12^1 in order to satisfy both these two features, 
short (£) and large C. Small-world networks have short-range 
connections between neighbors, as in regular networks, but 
they also have long-range connections similar to random net- 
works, without middle range ones. There are mainly two 
procedures to construct a small-world network: starting from 
same regular network, where each site is coupled to some 
number of first neighbors, one either rewires each regular 
connection with probability p (Watts-Strogatz model 12810 or 
adds a random connection to each node with probability p 
(Newman- Watts model l32lo . The second procedure is more 
appropriate for most purposes, since it avoids the possibility 
of generating disconnected clusters 13211 . 

Both random and small-world topologies produce typically 
networks where connections obey a Poisson distribution (see 
Tab. [J. However, there are real systems which are scale-free, 
i.e. where the connection distribution obeys a power-law. 

Scale-free networks were introduced by Barabasi and Al- 
bert l30ll using growing and preferential attachment: one starts 
with a small amount of nodes totally interconnected, and 
adds iteratively one node with m connections to the previ- 
ous nodes, chosen from a probability function proportional to 
their number of connections. With this construction one ob- 
tains analytically H |H |H a distribution of the connec- 
tions P(k) oc fc~ 7 , where 7 — > 3 as N — > 00, independently 
of the initial number of fully interconnected nodes and of m. 



It is also possible to generate scale-free networks, by either 
imposing a priori a power-law distribution of all connections 
randomly distributed, or by following a deterministic iterative 
rule for new nodes. The first procedure generates what is usu- 
ally called a generalized random graph, while the latter was 
recently referred as deterministic scale-free network l37ll . 

Deterministic scale-free networks, are hierarchical struc- 
tures composed by some succession of generations of nodes, 
i.e. the set of new nodes appearing simultaneously at a given 
iteration during the 'construction' of the network, whose con- 
nections follow a particular power-law distribution l37t I3H 
l39T.l40ll . being more easier to handle. The main difference be- 
tween random and deterministic scale-free networks is due to 
the local connectivity character of the latter: random construc- 
tions generate irregular long-range connections, while deter- 
ministic networks impose a succession of generation of new 
nodes which are, in some way, organized in 'space'. There- 
fore, deterministic networks are applied for instance in spin 
systems 13911 . and geographical and social networks I39t 14111 . 

After considering separately dynamical and topological 
complexity, the next logical step toward real network dynam- 
ics is to consider them together. One important question ad- 
dressed in this context is to know if full synchronization be- 
tween oscillators in such complex topologies would appear 
and under which conditions it is stable. By full synchroniza- 
tion we mean the convergence of the amplitudes of all oscil- 
lators to the same value, evolving coherently from then on. 
Therefore we call henceforth these fully synchronized states 
coherent states, to distinguish them from partially synchro- 
nized configurations, when several different clusters of nodes 
with the same amplitude are observed |16]. Synchroniza- 
tion and coherent behavior of oscillator networks with com- 
plex topologies have been studied for the random topo log y 
141 E3, HlH and small- world^ologyjES |3 El El, 
and also scale-free networks I4BI l50t IBH l52l 15311 . In random 
networks, it is already known [431 that with high coupling 
strengths it is possible to fully synchronize oscillators and 
the corresponding stability condition may be computed l43ll 
from the matrix of connections characterizing the network. 
In small-world networks, synchronizability is observed l47il 
only at the end of the small-world regime (high values of p), 
and recently 14611 it was found that heterogeneity in the cou- 
pling may destroy coherence. These findings for small-world 
networks are somehow controversial with the ones of other 
studies 1451 |48ll and other topological quantifiers have been 
proposed 14911 . In scale-free networks some recent studies 
indicate that synchronizability among oscillators depends on 
the average connectivity l52il and is robust to delayed flow of 
information l5(ill and to the removal of low-connected nodes 

In this manuscript we describe the general approach to 
study coherent states in any general complex network of os- 
cillators, and apply it to the particular case of a scale-free net- 
work of discrete-time oscillators, which is studied in great de- 
tail. We start in SectionllTlbv describing the stability analysis 
approach to model in Eq. Q and deduce the corresponding 
conditions for synchronizability. In Section|n]we apply this 
stability analysis procedure to the particular case of scale-free 
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networks, comparing our results with numerical simulations. 
The random scale-free case is treated in Section lTlI Al where 
we show that the threshold value of such a transition as a func- 
tion of coupling strength and outgoing connectivity obeys a 
power-law with an exponent that depends on the nonlinearity, 
while deterministic scale-free networks are studied in HUB I 
namel y a pseu do-fractal network l38ll and an Apollonian net- 
work I39H40I1 . Discussion and conclusions are given in Sec- 
tionEvl 

IL GENERAL APPROACH TO ANALYZE COHERENT 
STATES 

For all the network topologies described above, if one con- 
siders discrete-time oscillators, namely maps of the interval, 
the equation of evolution for their amplitudes reads 

&t+i = f{xt) - £&g{x t ), (1) 

where e is the coupling parameter, t labels time, x*t = 
(xt,i, ■ ■ ■ , iCt.Jv) with Xt,i representing the amplitude at time- 
step t of node i = 1, . . . , N, where N is the total number of 
nodes, / = (/(xi), . . . , f(x N )) and g = (g(xi), g(x N )) 
with / and g being real nonlinear functions, and G is the cou- 
pling (Laplacian) matrix, whose element Gij represents the 
relative strength with which node i is coupled to node j, and 

satisfies the conditions Ylj=i Gij = anc ^ Gu = 1 f° r a ^ 
i = 1, . . . , N. In general G is a non symmetric matrix. 

Usually, one choses g{x) = x when studying linear cou- 
pling, and g{x) = f(x) when studying nonlinear coupling. 
Here we consider the nonlinear case. Apart from this choice, 
all the information about dynamics is introduced in function 
f(x), while all the information about the coupling topology 
(regular, random, small-world or scale-free) and the coupling 
regime (either homogeneous or heterogeneous) is included in 
the coupling matrix G. 

From Eq. Q one easily sees that the coherent state xt.i = 
xt,2 = ■ • • = xt.N = Xt evolves in time according to the local 
map Xt+i = f{Xt). There are two ways to study these co- 
herent states: either by studying the stability of small pertur- 
bations of the coherent states or by making statistics over sig- 
nificant large samples of initial configurations, counting how 
many converge to a coherent state. Some attention to the pa- 
rameter ranges must be taken, since the basin of attraction of 
the coherent states may be bounded by regions of phase space 
where amplitudes diverge. In particular, for maps of the inter- 
val one has < e < 1 in order to guarantee convergence of 
any initial configuration. 

In this manuscript we will illustrate both analytical and nu- 
merical approaches for the particular case of scale-free net- 
works. To this end, we define the coupling matrix as Gu = 1 
and 

kf 



if node i is coupled to node j, with kj representing the number 
of nodes of node j and Kx is the set of labels of all neighbors 



of node i. If nodes i and j are not coupled Gy = 0. The pa- 
rameter a is a real number controlling the heterogeneity in the 
coupling: positive values of a enhance the coupling strength 
with sites having larger number of neighbors, while negative 
values favor sites having less neighbors. For a = the cou- 
pling between each site and its neighborhood is homogeneous. 

For local dynamics we choose the well-known quadratic 
map f(x) = 1 — ax 2 , where the free parameter a is restricted 
to the interval —0.25 < a < 2 and contains all possible dy- 
namical regimes from a fixed point (e.g. a = 0) to fully devel- 
oped chaotic orbits (e.g. a = 2). 

When determining the stability of coherent states, vari- 
ous criteria are possible. For instance, one could compute 
the maximum Lyapunov exponent and obtain the conditions 
where it is negative. However, maximum Lyapunov exponents 
do not indicate the existence of local instabilities in the syn- 
chronous state, which may pull the trajectories apart from the 
coherent manifold. 

The correct approach is based on the variational equation 
of Eq. Q proposed by Pecora and Carroll [:54|], which is valid 
for any network of identical oscillators in what concerns their 
local dynamics (quadratic map, Lorenz system, etc.) and their 
coupling regime (linear, nonlinear, etc). For the nonlinear 
coupling re gim e, the diagonal form of these variational equa- 
tions reads BlES El 

6+1., = exp(A(eXi))tt,i = [Df(X) - eXiDf(X)]Z t>i , 

(3) 

for coherent states xt.i = X, where A(eA,) is the Lyapunov 
exponent, Df(X) represents the identity matrix multiplied by 
the derivative of f(x) computed at x — X and Xi are eigen- 
values of the coupling matrix G. If G has zero-sum rows, 
i.e. Ylj=i Gij = OVi, and all its eigenvalues Ai < A2 < 
• • • < Xn are real and nonnegative, then Ai = corresponds 
to the mode parallel to the synchronization manifold and the 
largest Lyapunov exponent defines a master stability function 
15411 . The coherent state is stable whenever A(eAj) < for 

i = 2,. ..,7v insm. 

In our case, it is easy to check from Eq. that indeed 
G have zero-row sum, yielding Ai = 0. Moreover, all 
the eigenvalues of matrix G are real and non-negative, since 
det(G — AI) = det(G — AI) where G is a positive semidef- 
inite symmetric matrix, namely G = H 1/2 ]K 1 / 2 AK 1 / 2 II 1 / 2 
with A being the adjacency matrix of the network l52ll . and 
matrices H and K being the diagonal matrices with elements 

Ha = l/CEifjCj k t ) and K a = K respectively. 

From Eq. (|3|l and regarding the ordering of the eigenvalues 
Ai one easily concludes that the stability condition reads 

_ 1 - cxp (-A) 1 + exp (-A) _ 
£l = : < £ < = eu, (4) 

M Ajv 

where A is the Lyapunov exponent of the local single map. 
In particular there is a range of coupling strengths enabling 
synchronizability whenever Aat/ A2 < (1 + e~ A )/(l — e~ A ) 
holds. Therefore, by computing the eigenvalues of the Lapla- 
cian matrix G one is able to find the range of couplings for 
which coherent states are stable. This approach can be ap- 
plied for any system ruled by Eq. Q. 
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FIG. 1: Typical histograms of the standard mean square amplitude deviation a 2 showing the sharp transition to coherence as a function of 
(a) the coupling strength e with k = mo = 8, and (b) the outgoing connectivity k with e = 0.95. Values of N represent the fraction of 
configurations which converge to a coherent state from a total of 500 initial random configurations, after discarding transients of 10 4 time 
steps. Here a = 2, N = 1000, and a = 0. 
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FIG. 2: Boundary values el and eu in Eq. as a function of connectivity k (a) for a — 1.9 (dashed lines) and a = 2 (solid lines), and (b) 
for a = 1.5, 1.6, 1.7, 1.8, 1.9 and 2, whose inset is magnified in (c). The regions labeled with 'sync' are the ones where coherent solutions are 
observed, i.e. el < eu- Notice that in (c) the boundary for a = 1.8 is below the one for a = 1.7, contrary to other values (see text). 



III. SCALE-FREE NETWORKS OF COUPLED LOGISTIC 
MAPS: AN EXAMPLE 

For the particular case of scale-free networks, recent results 
14^1 show a transition to full synchronization for two partic- 



ular values of the nonlinearity a in the homogeneous regime 
(a = 0), when either the coupling strength or the number of 
outgoing connections are varied. However, as far as we know 
there is no detailed study showing how these coherent states 
depend on all the parameter models. Therefore, we present 
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in this Section detailed numerical results concerning synchro- 
nization in oscillator networks with scale-free topologies. Our 
purpose is to give a complete example of how to study coher- 
ent solutions in complex networks of oscillators, comparing 
both the stability analysis and the numerical approaches. 

The stability analysis is carried out just by computing the 
boundary values el and ejj in Eq. j§ as a function of the 
model parameters. The ranges of values where el < £u are 
the ones where coherent solutions appear. As stated above, 
for discrete oscillators, ruled by maps of the interval, the con- 
dition < £ < 1 must be added. 

Numerically, to detect coherent solutions from a given sam- 
ple of initial configurations, we compute the standard mean 
square deviation 14511 



N 

N 

i=l 



(x t ,i - x t y 



(5) 



where x t is the average amplitude at a given time step t. 
Whenever a 2 is zero within numerical precision, i.e. a 1 ~ 
10~ 30 , all the nodes are synchronized at the same amplitude. 

We divide our approach in two parts, the first one con- 
cerning random scale-free networks (Section BlI A> and a sec- 
ond one concerning deterministic scale-free networks (Sec- 
tion lTTTBt . 



A. Random scale-free networks 

In this Section, we use the algorithm of Barabasi and Al- 
bert 12, l30ll to construct the random scale-free network (see 
Section UJi, where at each node one places a chaotic logistic 
map. In a previous work l45ll a transition to coherence be- 
tween chaotic logistic maps was found for random scale-free 
networks, occurring for particularly high coupling strengths, 
typically of the order of e c ~ 0.9. Our simulations have 
shown that these transitions occur after discarding transients 
of ~ 10 4 time steps and they do not change significantly with 
the network size. Moreover, as shown in Fig.^ this transition 
to coherence is robust with respect to initial configurations, 
either by varying the coupling strength e (Fig-HJ) or the out- 
going connectivity k (Fig.[TJ>). In particular, above the thresh- 
old e c ~ 0.9, all initial configurations converge to a coherent 
state, indicating that in this parameter region the basin of at- 
traction of coherent states fills almost the entire phase space. 

From stability analysis, we find that in the fully chaotic 
regime (a = 2) the transition to coherence occurs for grad- 
ually smaller coupling strength if the connectivity k is in- 
creased. Figure [2^ shows the boundaries el and ejj as a 
function of k for a = 2 (solid lines) and a = 1.9 (dashed 
lines) with the same parameter values as in Fig. ^ As one 
sees, in both cases the lower boundary el decreases when 
k increases, while the upper boundary ejj increases beyond 
e = 1. Therefore, one expects that the region of synchroniz- 
ability increases for larger values of connectivity k. Figure 
|2k also shows clearly that for a — 2 the intersection between 
both boundaries, el = £u> occurs just above k = 7, which ex- 
plains why the transition to coherence in Fig.^ occurs at this 



value. For a = 1.9 this transition should occur near k = 5. 
Decreasing even more the nonlinearity coherent solutions are 
observed for even smaller connectivities and synchronizabil- 
ity regions increase, as shown in Fig. [2j>. To see this feature 
more clearly we magnify in Fig. |2j; the inset of Fig. |2j). As 
one sees, one exception occurs for a = 1.8, where the lower 
boundary is below the one for a = 1.7, due to the fact that for 
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FIG. 3: Transition to coherence as a function of connectivity k 
and coupling strength e. (a) Fraction Ncr = o of coherent states 
from 500 random initial configurations for a = 2. (b) Coher- 
ence transition curves in the (e, k) plane for (from bottom to top) 
a — 1.5, 1.6, 1.8, 1.7, 1.9 and a = 2, and (c) the same transition in 
a log-log plot, showing power-law dependence between connectivity 
and coupling strength for the transition curves, with an exponent /i 
which depends on the value of a (see inset). Here a = 0, L = 1000 
and we used transients of 10 4 time steps. By increasing the tran- 
sient size to ~ 10 6 one sees clearly that the transition to coherence 
is of first-order either (d) when varying the coupling strength e or (e) 
when varying the outgoing connectivity k. 
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FIG. 4: Boundary values el and Eu in Eq. J4j as a function of (a) connectivity k with a = 2 and a = (solid lines) a = 1 (dashed lines) and 
a = 2 (dotted lines), (b) nonlinearity a for fc = 8 and a = 0, 1 and 2, and (c) heterogeneity a with fc = 8 and o = 2 (solid lines) and a = 1.5 
(dashed lines). The inset in (c) emphasizes one small region where synchronizability is not observed, el > Eu (see text). Here N = 1000. 



a = 1.8 the Lyapunov exponent of the logistic map is smaller 
than the one for a = 1.7, as illustrated below in Fig.0] For all 
these values of a, the single uncoupled map shows chaotic or- 
bits. Moreover, for any other network size N, the same curves 
are obtained. 

These analytical predictions extracted from the stability 
condition in Eq. @ and shown in Fig.|2]are strongly corrob- 
orated with our numerical simulations as shown in Fig. [3] In 
Fig- EJi we plot the fraction F of initial configurations which 
converge to a coherent state for a = 2, while Fig.[3J> shows 
the threshold values, e c and k c , at the transition curves where 
the entire sample of initial configurations converge to a coher- 
ent state, for the same values of a as in Fig. |2j;. Here, one 
clearly sees that there is a clear and sharp transition to coher- 
ence. Interestingly, the curves in Fig.[5J> fit very well the ones 
in Fig. |2j;, which means that whenever the synchronizability 
condition el < £u holds, coherent states fill almost entirely 
the phase space. 

Moreover, as illustrated in Fig.|3J;, all curves obey a power- 
law, within our numerical precision, 

e c oc fc~". (6) 

For the six above values of a, the exponents are respectively 
H = 0.2345, 0.2354, 0.2353, 0.2231, 0.2023 and 0.1804: the 
exponent is almost constant below a ~ 1.7 and decreases 
above this value, as illustrated in the inset of Fig.[3J:. 

To determine the nature of the transition to coherence seen 
in Figs. [3^ andUJ), we plot in Figs. [JJi and[5£ the average 
standard deviation in the region where transition to coherence 
is observed, using much higher resolution. One clearly sees 
that the transition to coherence is of first-order, either when 
varying e or k. That the transitions are indeed of first order is 



easily recognized by the clear existence of hysteresis: when 
increasing either e or k the configuration eventually falls into 
a coherent state, no longer spontaneously desynchronizing, no 
matter how far the parameters are tuned back. 

All results till now, concern the case of homogeneous cou- 
pling (a = 0). Next, we study the case of heterogeneous 
coupling. Figure|4]shows the boundaries el and Eu in Eq. @ 
as a function of outgoing connectivity k, nonlinearity a and 
heterogeneity a, covering both the homogeneous and hetero- 
geneous regimes. Figure^ shows the two boundaries as a 
function of k for a = 2 and a = (solid lines), a = 1 (dashed 
lines) and a = 2 (dotted lines). As one sees for nonzero 
values of a the boundaries are no longer smooth curves, but 
instead they show fluctuations as k is increased, enlarging 
and shrinking alternately the region of synchronizability, la- 
beled as 'sync'. When varying a (Fig.^Ji) the boundaries are 
mainly controlled by the Lyapunov exponent of the local map 
(see Eq. @), where el (resp. eu) decreases (resp. increases) 
whenever a periodic window occurs l57ll . The fluctuations 
observed in Fig.0Ji are clearly seen in Fig.@J;, where the sta- 
bility boundaries are plotted as a function of a fixing k = 8 
and a = 2 (solid lines) and a = 1.5 (dashed lines). The fluc- 
tuations are much higher for a > 1 and for the fully chaotic 
regime both boundaries may even cross each other suppress- 
ing synchronizability (see inset of Fig. Moreover, the 
lower boundary el decreases till a ~ 0.5, then increases till 
a ~ 1 and decreases in average from there on. 

All these analytical results computed from Eq. @ and ma- 
trix G in Eq. Q are corroborated by our numerical simula- 
tions. In particular, the boundaries el and eu seen in Fig.[4j: 
are obtained also when plotting the contour of Fig. [5k where 
we plot the average standard deviation from a sample of 500 



FIG. 5: Transition to coherence as a function of the coupling strength e and heterogeneity a (a) when most connected nodes dominate the 
dynamics (a > 0) and (b) when the coupling to nodes with least neighbors is strengthened (a < 0). Here, we compute the average standard 
deviation from a sample of 500 initial configurations and fix a — 2, k = mo = 8 and N = 1000. 




FIG. 6: Transition to coherence when synchronization is imposed to all nodes having a number of neighbors (a) larger than a threshold kt, and 
(b) smaller than kt (see text). Here a = 2, k = 8, a = and N = 1000. 



initial configurations and vary the coupling strength and het- 
erogeneity for a = 2 and L = 1000. While Fig.[5k shows the 
numerical results for a > 0, i.e. in the case where nodes are 
more strongly coupled to the neighbors with higher connectiv- 
ities, Fig. [5J) shows the transition to coherence when a < 0. 
Here synchronizability is observed only for a > —0.15 and 
for very high coupling strengths e > 0.95. 



We end our study of coherent solutions in random scale-free 
networks by investigating briefly the role of hubs in the lattice. 
Instead of strengthening the coupling to the most connected 
nodes by increasing a > 0, we now fix a = and impose 
synchronization between all the nodes with more than a cer- 
tain threshold kt of neighbors and observe which fraction of 
the initial configurations converges to a coherent state. In this 
case the transition to coherence converges asymptotically to a 
limit of the coupling strength, as shown in Fig.|6^. The same 
occurs when synchronization is imposed to all nodes with less 
than kt neighbors, as shown in Fig.|6j5. 



B. Deterministic scale-free networks 

In the previous Section we focused on random scale-free 
networks, i.e. growing networks where new nodes are con- 
nected following probabilistic rules. In this Section we study 
deterministic scale-free networks l37ll38ll39Tl . using two dif- 
ferent deterministic topologies: the pseudo-fractal scale-free 
network introduced by Dorogovtsev et al I3q| and the Apol- 
lonian network introduced by Andrade et al 13911 and studied 
also in Ref. [4011. Both networks are illustrated in Fig.0 

The pseudo-fractal network of Dorogovtsev is obtained, 
starting from three interconnected nodes, and at each itera- 
tion each edge generates a new node, attached to its two ver- 
tices. Figure illustrates this network after three iterations, 
i.e. with three generations of nodes. The number of nodes N n 
and the number of connections V n increases with the number 
of generations as OF 



N n = |(3" + 1) 



Vn 



(7a) 
(7b) 



From Fig. [7^ one easily sees that this network has in- 
deed a scale-free topology, since the number of nodes 
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FIG. 7: Illustrations of two deterministic scale-free networks: (a) the 
pseudo-fractal network I4H1 . and (b) the Apollonian network 13^1 . 
Identical symbols label nodes belonging to the same generation n 
(see text), namely Q f° r n = 0, I for n = 1 and • for n = 2. 



with degree k = 2, 2 2 , . . . , 2"" 1 , 2™ and 2 n+1 is equal to 
3™, 3 n_1 , . . . , 3 2 , 3 and 3 respectively. In particular, the ex- 
ponent of this power-law distribution is 7 = 1 + In 3/ In 2. 
Moreover, the cluster coefficient of a node with degree k 
is C = 2/k, and the average path length is approximately 
(£) ~ 41niV n /(91n3). 

The Apollonian network is constructed in a different way: 
one starts with three interconnected nodes, defining a triangle; 
at n = one puts a new node at the center of the triangle 
and joins it to the three other nodes, thus defining three new 
smaller triangles; at iteration n = 1 one adds at the center 
of each of these three triangles a new node, connected to the 
three vertices of the triangle, defining nine new triangles and 
so on (see Fig. 05). The number of nodes and the number of 
connections are given respectively by 



N n = 

V n = 



i(3^ 
|(3"- 



5), 
1) • 



(8a) 
(8b) 



The distribution of connections obeys a power-law, since the 
number of nodes with degree k = 3, 3-2, 3-2 2 , . . . , 3-2 n_1 , 3- 
2" and 2 n+1 is equal to 3", 3 n ~\ 3 Tl " 2 , . . . , 3 2 , 3, 1 and 3 re- 
spectively, and the exponent 7 is the same as for the pseudo- 
fractal network. Moreover, a node with k neighbors has a 
cluster coefficient of C ~ 4/fc as reported in [40], converging 
on average to Coo = 0.828, and the average path length grows 
weaker than In Nr, (#%l . 

Despite that both networks have similar values for the topo- 
logical quantities, they are quite different from the geometri- 
cal point of view: the pseudo-fractal network has no metric, 
while the Apollonian network is embedded in Euclidean space 
and fills it densely as n — > 00, being particularly suited for de- 
scribing geographical situations 13911 . 

For stability analysis purposes (see Section ITU, the Lapla- 
cian matrix G of deterministic networks can be analytically 
determined from the adjacency matrix A = {a^ ■}, since they 
are related by 

G = I + AT, (9) 
where I is the identity matrix and the values of matrix T = 



{Tij} are defined by 



T 



dji 




ot 


Z^p=l a 3P 


Z^fe=l a pk 


Q 



(10) 



A simple way to write the adjacency matrix of the pseudo- 
fractal network is 



(11) 



N n xN n 



where N n is given by Eq. ( ITal . M T represents the transposed 
matrix of M and for each generation jt. = 1,2,... the matrix 
ML reads 



(12) 



2-3" 



X3" 



with 



A ... 
... A 



(13) 



1 vS»-l 



3™-!x3 



and whose starting form is 



^0 = 



1 1 

1 1 
1 1 



(14) 



3x3 



For the Apollonian network, the adjacency matrix is given 
by the same recurrence of Eq. (^J, but this time with 



111 
10 11 
110 1 
1110 



(15) 



and M n being a matrix with (3™ + 5) /2 rows and 3™ columns 
and having in each column only three non-zero elements. 

Figure [8] shows the eigenspectra of the Laplacian matri- 
ces for both the pseudo-fractal (Fig. |8^) and the Apollonian 
(Fig. IHJ?) networks, as a function of heterogeneity. As one 
sees for a = 2 (solid lines) synchronizability is observed only 
above a > 1.5, and in particular there is no synchronizabil- 
ity for the homogeneous coupling regime (a = 0). Figure 
|9ji shows the distribution of the average standard deviation 
over a sample of 500 initial configurations, from which one 
clearly sees that there are no coherent solutions. Here the 
standard mean deviation is characterized by some large value 
which is almost constant beyond the weak coupling regime 
(e > 0.2). In the weak coupling regime (e < 0.2) the stan- 
dard mean square deviation is even larger, since the coupling 
is not strong enough to compensate the highly chaotic local 
dynamics (a = 2). 

From Fig. [8] one also sees that, for the pseudo-fractal and 
a > 1.5, the upper threshold e\j increases monotonically with 
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FIG. 8: Boundary el and eu for synchronizability as a function of heterogeneity a for (a) the pseudo-fractal network and (b) the Apollonian 
network, with a = 2 (solid lines) and a — 1.5 (dashed lines). For each network we use 6 generations of nodes (see text). These eigenspectra 
are the same for any number of generations. 




FIG. 9: (a) Typical histogram of the standard mean square deviation a 2 for the pseudo-fractal network as a function of the coupling strength 
e, with a — 2 and a — 0. A similar result is obtained for the Apollonian network, (b) Histogram of the standard mean square deviation a as 
a function of nonlinearity a and coupling strength e, for deterministic scale-free networks with a — 0. The mean square deviation is averaged 
over a sample of 500 initial configurations and during 100 time steps, after discarding transients of 10 4 time steps. 



the heterogeneity, while for the Apollonian network the upper 
threshold decreases. This particular difference between both 
networks should be due to their geometrical differences, in 
particular the fact that Apollonian networks are embedded in 
Euclidean space could explain in some way that stronger dom- 
inance in the coupling to the most connected nodes destroys 
coherence. 

Choosing other values of a for which local dynamics is 
chaotic, one finds the same values of el and sjj as functions 
of a only shifted: sl gets smaller, while eu increases. Fig- 



ure [8] illustrates this for the particular case of a ~ 1.5. De- 
creasing even further the nonlinearity below the accumulation 
point a = 1.411 . . . synchronizability is attained for any pos- 
itive value of a, whenever the coupling strength is sufficiently 
strong. 

Moreover there is a complicated dependence of the aver- 
age standard deviation on the coupling strength and nonlin- 
earity. As shown in Fig. |9p for deterministic scale-free net- 
works one finds two main regions in the (a, e) plane: (I) a 
region where the standard mean square deviation is large and 
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FIG. 10: Inducing transition to coherence by varying the heterogeneity a (see Eq. 0) in scale-free networks, (a) pseudo-fractal network and 
(b) Apollonian network. For strong heterogeneity coherence appears beyond a relatively high coupling strength, and disappears again for very 
large couplings (see text). For each network, we use £ = 6 generations of nodes and fix a = 2. (c) and (d) show high-resolution plots of a 2 as 
a function of e for a = 2, emphasizing the first-order phase transition to coherence. 




FIG. 1 1 : Transitions to coherence in deterministic scale-free networks, when synchronizing the first g generations of nodes out of I generations 
(see text), (a) pseudo-fractal network and (b) Apollonian network. The collective dynamical behavior is quite insensitive to hubs (see text). 
Insets show that transitions to coherence are of first-order. For each network, we use £ = 9 generations of nodes and a — 2 fixed. 



varies smoothly with the parameters and (II) a region where 
the mean square deviation is smaller but has larger fluctua- 
tions. This second region, observed for a < 1.7, is some- 
how surprising, since irregular variations of the standard mean 
square deviation occur for low nonlinearity and high coupling 
strengths, precisely where one would expect the most regular 
behavior of the node dynamics. 



As for the heterogeneous coupling regime (a ^ 0), Fig. 1101 
illustrates the transition to coherence by varying the hetero- 
geneity a for the pseudo-fractal (Fig.llOh) and the Apollonian 
network (Fig. 110b). For both networks, one sees that coher- 
ence sets in for a > 1.5, and the contour of the histograms 
marking the transition to coherence fits well the regions in 
Fig. |3 labeled as 'sync'. Moreover, from Figs. 110b and 1 1 OB 
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one observes that all these transitions to coherence are of first- 
order. 

Finally, we study the role of hubs in deterministic scale- 
free networks, as we did in the previous Section for random 
networks. To this end, we impose synchronization among 
g = 1 , . . . , I generations, with £ being the total number of 
generations, and observed in what conditions coherent states 
are observed. In the pseudo-fractal network the first genera- 
tion has Ni — No = 3 nodes, the second one has N2—N1 = 9 
nodes, and the nth generation has N n — N n _i = 3™ nodes. 
In the Apollonian network the number of nodes appearing at 
each generation is precisely the same. 

Figure fTTI shows the standard mean square deviation as a 
function of coupling strength for pseudo-fractal (~Fig.ll lh) and 
Apollonian networks (Fig. II lb). In each case we choose the 
fully chaotic map (a = 2) and impose synchronization among 
the nodes of the first g generations by setting them to be their 
mean amplitude at each time-step. In both cases, one sees that 
the standard mean square deviation remains large when syn- 
chronization is imposed to all g < 1 — 2 generations. Coherent 
solutions are only observed for g = £ — 2 and g = I — 1, be- 
yond a coupling threshold which is smaller for the latter case. 
Surprisingly, for g = £ — 1 the transition to coherence occurs 
for the same coupling strength in both networks. This may 
be due to the fact that, the fraction N g /N e = 3 9 ~ £ = 1/3 of 
nodes on which one imposes synchronization is the same for 
both networks and is high enough to suppress the influence of 
local connectivities. 

For g — £ — 2 the pseudo-fractal network shows coherence 
only above very high coupling strengths, near e <~ 1, while 
for Apollonian networks the threshold is much lower. This 
difference in the coupling strength threshold is due to the fact 
that here the fraction of nodes N g /Ni = 1/9 to which syn- 
chronization is imposed is small enough to not suppress the 
influence of local connectivities. Therefore, since the hubs in 
the pseudo-fractal network are less connected than the hubs 
in Apollonian networks, one needs higher coupling strength 
to observe coherence. For any higher value £ of generations 
the same results are obtained, since one has for the quotient 
of the number of nodes between two successive generations 
N n /N n -i — > 3 as n increases. 

As a general remark, one observes from Fig. ^2 that one 
needs to synchronize a rather high fraction of nodes (> 1/9) 
to induce coherence. Therefore, it seems that, dynamical col- 
lective behavior on scale-free networks is quite insensitive to 
hubs. As shown in the insets of Figs. II lh andll lb. the transi- 
tion to coherence is also of first-order. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we studied fully synchronized solutions for 
three scale-free network topologies. The main conclusion is 
the following: in random scale-free networks synchronization 
of chaotic maps not only depends on the coupling strength but 
is mainly controlled by the outgoing connectivity k, which 
is a measure of cohesion in the networks. Because of that, 
one finds coherent solutions in random scale-free networks of 



fully chaotic logistic maps (a = 2) with outgoing connectiv- 
ity k = 8 and homogeneous coupling, but not in determinis- 
tic scale-free networks, since they have rather small effective 
outgoing connectivity, namely k = 2 for the pseudo-fractal 
network and k = 3 for the Apollonian network. Therefore, 
although the exponent 7 of connection distributions in scale- 
free networks does not depend on the outgoing connectivity 
we have shown that, in general, synchronization of chaotic 
maps in such coupling topologies is quite sensitive to it. 

Our results were obtained both numerically, from histogram 
of significantly large samples of initial configurations with a 
criterion for full synchronization based on the mean standard 
deviation of amplitudes, Eq. Q, and analytically from the 
eigenvalue spectra of the diagonalized variational equations 
computed at the coherent states, Eq. 0. 

In particular, for random scale-free networks, the threshold 
values of the coupling strength obey a power-law, Eq. ©, as 
function of the outgoing connectivity. The exponent of this 
power-law depends on the nonlinearity a of the chaotic map, 
being almost constant below a c ~ 1.7 and decreasing linearly 
above it. Interestingly this value of a c is in the vicinity of the 
bifurcation of the quadratic map where the period-3 window 
appears, and coincides with the appearance of other nontriv- 
ial behaviors in coupled map lattices with regular topologies, 
namely in the velocity distribution of traveling wave solutions 
0. 

For deterministic scale-free networks with homogeneous 
coupling, the same value a c indicates the threshold above 
which no coherent solutions are observed, independently of 
the coupling strength. Above a c , coherence is observed only 
for heterogeneous coupling, namely for a > 1.5. However, 
for this range of values, we have also shown that coherence 
is also absent either for very small or for very large coupling 
strengths, due to spatial instabilities. Another particularly in- 
teresting result that still needs to be explained is that, for Apol- 
lonian networks, the coupling threshold beyond which coher- 
ence disappears gets smaller when the heterogeneity is further 
increased. This point is not observed for the pseudo-fractal 
network and may be due to the geometrical differences be- 
tween both deterministic networks. 

As a general property, we have shown that all transitions to 
coherence are of first-order, indicating a similarity with other 
complex networks l2*jl . Furthermore, all results are robust 
not only against changes of the initial configurations of node 
amplitude but also, in random scale-free networks, against 
changes of the connection network. We also presented results 
indicating that in scale-free networks hubs play apparently no 
fundamental role in the dynamical collective behavior. 
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